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Abstract 

We study dilepton production from a thermally equilibrating quark-gluon 
plasma expected to be formed in ultra relativistic heavy-ion collisions. The 
pre-equilibrium dynamics of quark-gluon plasma is studied within the color 
flux-tube model by solving non-abelian relativistic transport equations. This 
dilepton rate crucially depends on the collision time of the plasma. The 
results are compared with the Drell-Yan productions. We suggest that a 
measurement of this rate at RHIC and LHC will determine the initial field 
energy density. 
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One of the primary aims of the ultra relativistic heavy-ion collision experiment is to detect 
a new phase of matter, viz quark-gluon-plasma(QGP). It is known from lattice quantum 
chromodynamics that hadronic matter undergoes a transition to this quark-gluon plasma 
phase at very high temperatures (~ 200MeV) and densities (~ 2GeV/ fm 3 ) Q. While 
such a phase did exist in the early universe, after ~ 10~ 4 seconds of the big bang, it is 
interesting if we can recreate the early universe and experimentally verify this QCD phase 
transition in the laboratory i.e. in high energy heavy-ion collisions. In the near future, the 
relativistic heavy ion collider(RHIC) at BNL (Au-Au collisions at = 200GeV) and large 
Hadron Collider(LHC) at CERN (Pb-Pb collisions at ~ 6ATeV) @,|| will provide the 
best oppurtunity to study such a phase transition. As this plasma exist for a very short 
time (~ 10 fm) and formed over a very small volume (~ 100/m 3 ), the direct detection of 
this phase is not possible in these experiments. The proposed signatures for its detection are 
therefore indirect, and the prominent ones are : 1) J/ip supression @], 2) electromagnetic 
probes such as dilepton and direct photon production |5|,|| , and 3) strangeness enhancement 

It has been suggested by Matsui and Satz that J / ip suppression is a good probe (3| for the 
detection of an equilibrated quark-gluon plasma. Their calculation is based on lattice QCD 
which assumes equilibration of the plasma. However, in heavy-ion collisions such as at RHIC 
and LHC the assumption of such an equilibration of quark-gluon plasma is only suggestive. 
One does not know exactly when the plasma equilibrates and there is no concrete calculation 
on this aspects. This sets an uncertainty in the calculation of J /if) survival probability which 
is determined from screening. In this regard, there are calculations of J /if) suppressions in 
the equilibrating quark-gluon plasma HH using short- distance QCD, but the uncertainties 
in the J /if) formation time and quark-gluon plasma equilibration time make it difficult to 
compare the results with the experiments (see Ref. f| for the details). The suppression of 
J /if) infact, been observed in reactions where there is no quark-gluon plasma phase, such as in 
p- A collisions and in collisions of light nucleus . This suppression is well explained by the 
nuclear absorption of J /if) [pjj . Although, recently, NA50 collaboration |12|] reports an excess 



in the supression of J/ip, it is still not clear if an equilibrated quark-gluon plasma has formed 
in this collisions. There are proposals that the data are explained by a deconfmed partonic 
medium |y|, or by a medium with high density flijl ; but there are also other calculations 



which explain the data without assuming any quark-gluon plasma phase |15|,|16|. As far as 
the J/ip supression in nucleus-nucleus collision is concerned, many aspects of it has to be 
studied in greater detail before unambigious conclusions can be drawn about the existence of 
quark-gluon plasma. It thus appears that there is a need to study more than one signature 
if one has to detect QGP. With that in mind we investigate another signature, i.e. dilepton 
emission. 

Dileptons and single photons have long been proposed as useful probes of the plasma || , 
as once produced, they hardly interact with the strong matter and thus carry the details of 
the circumstances of their production. Being electromagnetic in nature, they do not suffer 
from the final state ineractions (interactions with hadrons) and hence keep the memory of 
their formation surroundings (such as the information about the temperature and number 
density). The major processes for dilepton production in ultrarelativistic heavy- ion colli- 
sions (URHIC) are, (i) hard parton scattering, (ii) electromegnetic decay of hadrons and 
(iii) production from partons present in QGP (thermal production). The hard parton scat- 
terings produce high p t lepton pairs which are calculated from pQCD, knowing the structure 
functions of partons inside nucleus. On the other hand, electromagnetic decay of produced 
hadrons is the main source of dilepton and photon production in heavy-ion collisions which 
obscure the signal of interest (thermal production), which are produced from QGP. The ther- 
mal emission of high invariant mass dilepton is calculated in the same way as that of hard 
parton scatterings by using the thermal distributions of partons inside QGP. It is expected 
that in high energy heavy-ion collisions such as at RHIC and LHC, the thermal production 
of dileptons will be more than that from other processes ||17|| . On the experimental side, 
an enhancement of dilepton yield is observed in central 200A GeV S+Au [fy|, S+W fll9f , 
and 160A GeV Pb+Au |2(J collisions. However, these data are partially explained by the 
conventional mechanism of binary hadron collisions, e.g., by tttt — > e + e~ processes [f21|-|23| 



and by the contribution of the collective bremsstrahlung mechanism In ultrarelativistic 



heavy-ion collisions such as at RHIC and LHC situation might be different and one may 
expect an enhancement of the dilepton production which can be due to the formation of 
quark-gluon plasma. 

For a thermally equilibrated quark-gluon plasma, the distribution of dilepton is given by 



where T c is the freeze-out temperature at which one stops the hydrodynamic evolution of 
the quark-gluon plasma, or perhaps the critical temperature if there is a phase transition. 
Here M is the invariant mass of the dilepton pair, R is the size of the nucleus and a is the 
anihilation cross section for quark-antiquark into a dilepton pair. The function H is written 
in terms of Bessel function K\ as, H(z) = z 2 (z 2 + 8)Ko(z) + Az(z 2 + 4)Ki(z). Assuming 
a fast thermalisation and high initial temperature, enhancement in the dilepton rate was 
found by several authors |17|]. However, it can be seen from the above equation that there is 
always an uncertainty in the dilepton rate because of the assumption of initial temperature 
(To) and initial time (tq), at which quark-gluon plasma thermalise. The dilepton rate in 
any realistic calculation depends crucially on these initial conditions which determines the 
hydrodynamic evolution. For this reason a detailed calculation of the plasma evolution in 
different stages of ultrarelativistic heavy-ion collisions is necessary. The various stages by 
which the evolution of quark-gluon plasma is described in URHIC are, i) pre-equilibrium, 
ii) equilibrium, where one actually studies the equilibrated quark-gluon plasma, and iii) 
hadronisation. The pre-equilibrium stage of the collision which leads to thermal and then 
chemical equilibrium has a crucial role to play in the equilibration of the plasma and also on 
the calculation of different signatures. ^From this point of view it is necessary to study what 
happens to the dilepton production from different stages of QGP, rather than estimating it 
in an equilibrated quark-gluon plasma (Eqn.-l ). In this paper we study the dilepton rate 
for a thermally equilibrating quark-gluon plasma in ultra relativistic heavy-ion collisions. 
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The rate of equilibration of quark-gluon plasma in URHIC is different for diiferent models. 
One of the relevant model that describes the production and the equilibration of QGP in 
URHIC is the color flux-tube model p5-27[. This model is a generalisation of the familiar 
Lund string model widely used for e + e~ and p — p collisions [pSfl . Within this model, two 
nuclei that undergo a central collision at ultra high energies are highly lorentz-contracted 
as thin plates. When these two highly lorentz-contracted nuclei pass through each other 
they acquire a nonzero color charge (< Q >= 0, < Q 2 >^ 0), by a random exchange of 
soft gluons. The nuclei which act as color capacitor plates produce a chromo-electric field 



between them [29,30] . This strong electric field creates qq and gluon pairs via the Schwinger 



mechanism which enforces the instability of the vacuum in the presence of an external 
field. The partons so produced, collide with each other and also get accelerated by the 
background field. In the case at hand, the color degree has a central dynamical role in the 
evolution of the plasma. In our recent studies we have incorporated this dynamics in the 
plasma evolution in heavy-ion collisions [|3^,[£|. The color charge(Q a ) which is a vector in 



the color space obeys Wong's equation |34 

dQ a 



dr 



f abc u^Q b A c ^. 



(2) 



where A a>1 is the gauge potential, and f abc is the structure constant of the gauge group. 
This equation describes the precission of the color charge in external chromo field. The 
non-abelian version of the Lorenz force equation 

dp^ 



dr 



Q a F a " u u u , 



(3) 



describes the acceleration of the color particles by the background chromo field F a ^ u . Due to 
these dynamical nature of the color charge, the usual classical phase space of coordinate and 
momenta is extended to include color. The single particle distribution function f(x,p, Q) of 
quark and gluon is then defined in the compact phase space of dimension 14 in SU(3). In 
this extended phase space a typical relativistic transport equation is written as |]35| 



V ^ + Q a F;yd; + r hc Q a Ayd c Q f(x, P , q) = c(x, P , q) + s(x, P , q). 



(4) 



One needs to solve this transport equation to study the non-equilibrium dynamics of the 
quark-gluon plasma expected to be formed in ultra relativistic heavy-ion collisions. The 
first term in the above non-abelian relativistic transport equation corresponds to the usual 
convective flow (free streaming expansion), the second term is the non-Abelian version of 
the Lorentz force term and the last term corresponds to the precession of the color charge 
as described by Wong's equation. S on the right hand side of equation (4) correspond to 
the source term for parton production from background chromoelectric field via Schwinger 
non-perturbative mechanism. C corresponds to the collision term. In general we have to 
write separate equations for quarks, antiquarks and gluons since they belong to different 
representations of the gauge group. For anti-quarks the distribution function f(x,p,Q) 
obeys a simillar equation, with Q a replaced by — Q a [i.e. the second term in the above 
equation changes sign). These equations, which are lorentz and gauge invariant ||35|| , are 
closed with the Yang-Mills equation, 



(D^) a (x) = f{x) = g J p v Q a [f q (x^Q) - f q (x,p,Q) + f g (x,p,Q)]dPdQ. (5) 

to study the production and equilibration of quark-gluon plasma. 

According to Bjorken's proposal [|36| the distribution function and other physical quan- 



tities are written in terms of the boost invariant parameters r(= \Jt 2 — z 2 ), p t and 
£(= 7] — y). Here rj is the space time rapidity (tanhrj = z/t) and y is the momentum 
rapidity (tanhy = p z /p°). In our calculation we have employed a collision term in the 
relaxation time approximation: 

c = -*■"»(/-/->, (6) 

where f eq is the local eqilibrium distribution function with explicit color dependences. This 
is given by [|37| : 



f eq = • (7) 



'" " exp((pM - Q a ■ A» a )ujT(r)) ± 1 ' 



The +(-) sign in the above expression is for quarks (gluons). The source term S for qq 
and gluon pair production is obtained by the Schwinger mechanism of particle production 



53fl . In the process where field and particles are present, we use the energy momentum 
conservation equation: 

9fj,T matter + d^Tj ield = 0, (8) 

with TZtter = Jp^(2f q + 2f q + f g )dTdQ and T^ eld = diag(£ 2 /2, E 2 /2, E 2 /2, -E 2 /2). Here 
dT = d 3 p/ (2ir) 3 p = p t dp t d!; / (2n) 2 , and dQ is the integral in the color space. The factor 2 in 
the above expression is for two flavors of massless quark. Equation (4) along with equation 
(8) is solved numerically by a double self consistent method |38j to study the equilibration of 
quark-gluon plasma in URHIC. In the earlier paper |33| we have studied the bulk properties 
and the equilibration of quark-gluon plasma (which occurs arond 1 fm). Here we calculate 
the dilepton production rate using these distribution functions of partons /(t, £,p t , Q). 
The dominent process which produces a dilepton pair is 

q + q^1*^l + + l-, (9) 

where 7* is the intermediate virtual photon. The dilepton emission rate dN for such a 
process in a space time volume d A x is 

dN 1 



■ J d 3 Pl d 3 p 2 dQf(x, Pl , Q)f(x, pa, Q)v rel a(M 2 )5\P - Pl - p 2 ). (10) 



d 4 xd 4 P (2tt) 6 

Here P M is the four momentum of the lepton pair, Pt is the transverse momentum, Mt( 



M 2 + Pt ) is its transverse mass and M is the invariant mass (M = P^P^). v re i(= 
M 2 /2EiE 2 ) is the relative flux velocity of quark and antiquark pair in the above process. 
The dilepton production cross section a(M 2 ) for the above reaction is 



.„ Ana r_ 2mf nr , 4mf n „ ,,,, 
a{M)= 3NP [1 + l^ ][1 -M^ ]F <' (U) 

with F q = N 2 jj J2f e } an d mi is the mass of the lepton. Here, N is the color averaging factor 
(which corresponds to the volume of the color space), N s is the spin degeneracy (iV s = 2s + 1) 
and ef is the fractional charge of the flavour. For dilepton pair of large invariant mass the 
emission rate in the midrapidity region (Y = 0) is given by: 



dN 5R 2 a 2 



dM T dYdP% 72vr 7 



drrW{hJ 2 ). (12) 



where 



urtf n f + °° a f + °° ac f P+ a f An Ptif( T iPti,£uQ)f(r,pt2,&,Q) n ^ 
W( fl , h ) = drj d£i I d Pa J dQ ^ ,, ( (13) 



\p 2 n P 2 - \p tl M t ch(rj - 6) - |M 2 ] 2 ]V2 



with p t2 = <JM T - 2M T p n ch(rj - £i) + Pa, sh£ 2 = ^(M T shr) - p tl sh£i) and p± = 
^M 2 [M T ch(r] — ^x) =)= Pt] _1 - in equation (12) we have used d 4 x = 7iR 2 dTrdr] (where R 
is the radius of the nucleus). a{= 1/137) is the coupling constant of the electromagnetic 
interaction. 

It can be mentioned here that the thermal dilepton rate does not depend on M and Pt 
separately, but only on M T [J39||. However this scaling is violated for an equilibrating 



plasma, which is seen in equation (12). In this case the dilepton rate depends on both 
and Pt- This would also be the case if one takes the transverse expansion of the plasma 



into account 40 



We will now present our results of the dilepton rate using the above non-equilibrium dis- 
tribution function, f(r, £,pt, Q), of quarks and antiquarks. In our calculation we have taken 
R = 7fm and the initial field energy density eo(= 1/2Pq) to be equal to 300 GeV//m 3 (see 
Ref. [Q). We have considered two different cases corresponding to two different relaxation 
times, r c = 5/m and r c = 0.2/m. The relaxation time r c = 5fm corresponds to collisionless 
limit and r c = 0.2 fm corresponds to a more realistic limit of equilibration ||41|| . 

We have calculated the dilepton rate as a function of for P? = and 1 GeV respec- 
tively. In Fig-1 we have presented the results for r c = 0.2/m. As can be seen from the figure, 
the dilepton yield becomes smaller for higher values of transverse momenta. In Fig-2 we 
have presented our results for r c = 5fm which corresponds to the collisionless limit. In this 
case the rate at higher Mt is found to be larger than that at r c = 0.2/m. This is because 
the average energy per parton is higher at r c = 5fm than at r c = 0.2/m as observed earlier 
f33| . The average energy per parton is around 4 GeV for r c = 5 fm and is around 2 GeV for 
t c = 0.2 fm. In the absence of any collision the partons come with higher energies than the 



partons with collision. From this point of view it is crucial to determine the collision time 
accurately. We will mention here that a determination of the collision time r c is possible 
by a triple self consistent numerical method [42| instead of a double self consistent method 



(which we have employed here to solve the transport equations). 

In order to get a feeling, how crucial the dilepton production from an equlibrating QGP 
is, we have compared our results with the Drell-Yan productions. The Drell-Yan results are 
taken from the Ref. [ji3|| . These spectra are calculated at pp center-of-mass energy a/s = 
200 GeV scaled to UU central collisions assuming a simple factorisation of the nuclear mass 
dependences. It can be seen in Fig-3 that the dilepton rate from the pre-equilibrium stage 
is larger than the Drell-Yan production for very small transverse momentum (P? — OGeV) . 
This is true when Mt < 2GeV. For > 2GeV, the Drell-Yan production dominates over 
the pre-equilibrium dilepton production. For large transverse momentum (Pt = 1 GeV) 
of dilepton pair, the production from the pre-equilibrium stage is smaller than the Drell- 
Yan production in the whole range of dilepton transverse mass. This has been shown in 



Fig-4. This result is contrary to the earlier findings |44],f|§] where it is shown that dilepton 
production from the pre-equilibrium stage dominates over the Drell-Yan production. The 
partons formed are not too hard in the color flux-tube model, and one expects to have such a 
low rate in the dilepton spectra than the Drell-Yan emissions, which are produced from the 
primary hard scattering of partons. However, the comparisons are not very strict because the 
choice of initial field energy density eo (= 1/2 Eq, Eq being the initial chromoelectric field) 
is arbitary and there is no way of determing this. This is because one does not know how 
many soft gluons are exchanged (which deterimines the strength of the initial chromoelectric 
field) when two nuclei cross each other in ultra relativistic heavy-ion collisions. This initial 
condition can be determined from the measurement of some experimental signatures, such 
as dilepton production. However, the results we have found here do not seem unnatural. 
This is because, if one dentines the energy deposited at RHIC |3(J (when two nuclei are 
in maximum overlap) with the initial field energy, the initial field energy density becomes 
~ 250 GeV/fm 3 [[|T|. In any case, we hope to extract the inital field energy density from 



the measurement of dilepton spectra at RHIC and LHC (this is done only after a careful 
substraction of the dilepton rates from other processes in different stages of quark-gluon 
plasma). After the determination of this initial field energy density, other bulk properties 
of the plasma can be determined accurately. 

Summarizing the paper, we have calculated the dilepton spectra in ultra relativistic 
heavy-ion collisions within color flux-tube model, with non-abelian features explicitly in- 
corporated. The production is larger in the collisionless limit of the plasma. After one 



determines the collision time accurately by a triple self consistent procedure [42 



an ex- 



perimental measurement of the dilepton spectra at RHIC and LHC will shed light on the 
determination of initial field energy density which we have assumed in our model. Only 
after this, the predictons of all the bulk properties of the plasma such as temperature and 
number density, will be determined accurately. 
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Figure captions 

FIG. 1. Dilepton rate from pre-equilibrium stage, as a function of for r c =0.2 fm. 
FIG. 1. Dilepton rate from pre-equilibrium stage, as a function of for r c =5.0 fm 
(collisionless limit). 

FIG. 3. Drell-Yan and pre-equilibrium dilepton rate as a function of for r c =0.2 fm 
(P T = GeV). 

FIG. 4. Drell-Yan and pre-equilibrium dilepton rate as a function of for r c =0.2 fm 
(P T = 1 GeV). 
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